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ORBITAL EVOLUTION OF COMPACT WHITE DWARF BINARIES 
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ABSTRACT 

The new-found prevalence of extremely low mass (ELM, Mh c < 0.2 M ) helium white dwarfs 
(WDs) in tight binaries with more massive WDs has raised our interest in understanding the nature 
of their mass transfer. Possessing small (M env ~ 1O~ 3 M ) but thick hydrogen envelopes, these 
objects have larger radii than cold WDs and so initiate mass transfer of H-rich material at orbital 
periods of 6-10 minutes. Building on the original work of D'Antona et al., we confirm the 10 6 yr 
period of continued inspiral with mass transfer of H-rich matter and highlight that the inspiraling 
direct-impact double WD binary HM Cancri likely has an ELM WD donor. The ELM WDs have less 
of a radius expansion under mass loss, thus enabling a larger range of donor masses that can stably 
transfer matter and become a He mass transferring AM CVn binary. Even once in the long-lived AM 
CVn mass transferring stage, these He WDs have larger radii due to their higher entropy from the 
prolonged H burning stage. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: Type la — white 
dwarfs 



1. INTRODUCTION 

Helium core white dwarfs (WDs) are made from 
< 2.0 Mq stars when stellar evolution is truncated 
before the He core reaches the Mh c ~ 0.48 M needed 
for the helium core flash. One formation mechanism is 
significant mass loss due to stellar winds o n the red giant 
branc h (RGB) that strips the H envelope (iD'Cruz et al 



1996), typically leading to M H « = 0.4 — 0.48 (iHanseii 



20051: [Kalirai et all l2007t iKilic. Stanek. k Pinsonneault 
2007bf K Another mechanism is a common envelope 



induced by binary interacti ons ()Iben k Livid 119931 : 
iMarsh. Dhillon. k Duckl 119951 ) , making extremely low- 
mass (ELM) He WDs (M He < 0.20 Af© or so) when 
the interaction occurs at the base of t he RGB (see 
Ivan Kerkwiik. Bergeron, k Kurkarnilll996l) . These ELM 
He WDs were first seen as comp anions to millisecond 
pulsars (e.g., Ba ssa et"al"1 120061) or in high proper 
motio n catalogs (jKawka et al . 2006; Kawk a k Vennesl 
120091 ) . but the advent of the Sloan Digital Sky 
Survey (jEisenstein et al.l I2006D and other surveys re- 
vealed many additio n al ELM WDs (IKilic et all 12 007b: 



Badenes et ail 120091: iMullallv et all 120091: IKilic et al 
2010bt IMarsh et al.l 120111: Kulkarni k van Kerkwiik 



20ToV ISteinfadt et alJ I2010bt IKilic et al.l 1201 Ida 



Parsons et all 120111: IMarsh! 120111: IBrown et al.l 12011b : 
Kilic et a.1.1 12011bt iVennes et al.l 120111 : IKilic et al.l 12012 . 
and Figure [T]). 

ELM WDs were predicted to possess stably burn- 
ing H envelopes {M m v ~ 10" 3 - 1Q- 2 M ) that keep 
them bright for Gyrs (|Serenelli et all l2002t iPanei et all 
120071 ) , and this has cer t ainly aided the re cent detec- 
tions (jKilic et all I2011ah IBrown et all 120101 ) . Though 
identified by their location in \ogg — T e g space, few 
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systems have actually had their radii m easured with 
any precision. Steinfadt et al's (|2010bl ) discovery of 
the eclipsing double W D binary system NLTT 11748 
(jKawka k Vennesl 120091 ) allowed for the first geomet- 
ric measurement of the radius of a ELM WD, finding 
R fa 0.04 R Q for the 0.15 Mq He WD, consistent with 
the prese nce of a thick stably- burning hydro g en envelope 
(also see iKawka. Vennes. k Vaccard" 120101: IKilic et al ' 



2010a ). Additional eclip sing systems (|Parsons et al 
2l CTiBrown"et~aTl 12011b! ) have led to even more con 



straints, although for some of the more compact systems 
it is not clear if the radius is truly the equilibrium radius 
of the WD or if it has been tidally distorted. 

Multiple common envelope phases are possible in the 
formation scenario for ELM WDs, leading to ELM 
WDs in binaries with more massive WDs. Many 
of these will come into contact within lOGyr (Fig- 
ure [T]). Indeed, a large number of the known dou- 
ble WD binaries contain ELM WDs (M He < 0.20 M , 
circled points) with large (> 0.0 3 R@} radii indica- 
tive of a s table H burning shel l (IKawka et all 120061: 
Kilic et all l2007al: iBadenes et all 120091: IMullallv et all 
2009t IKilic et alj|2010bt iKulkarni k van Kerkwiikll2010l: 
Marsh et al.ll2011t IBrown et alJl2011bHKilic et alll2011c 



Vennes et alll201lD . As noted bv lD'Antona et all (|2006D 



since the time to burn the H envelope can easily exceed 
the time to reach contact, many of these ELM WDs will 
come into contact with the remaining H envelope. This 
raises the possibility for many new phenomena that we 
begin to explore here. 

The large radii of the ELM WDs means that Roche 
lobe overflow (RLO) occurs at larger orbital periods than 
otherwise expected, so we start in §2 by examining the 
behavior of the radius of the ELM WD as its H envelope 
is transferred. We follow in §3 by outlining the basics 
of mass transfer and highlighting some of the new pos- 
sibilities when ELM WDs are donors. There is a more 
distinct possibility for thermally stable burning of the 
accreted H and He on the accreting WD and the initial 
contraction of the ELM WD to mass loss allows for more 
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Fig. 1. — The population of double WDs (triangles and squares, 
depending on whether one or two members of the binary have ra- 
dial velocity measurements), sdB/WDs (filled diamonds), and pul- 
sar /WDs (filled stars) with P or b < day. Binaries to the left of the 
dashed lines will merge in less than 10 Gyr, 100 Myr, or 1 Myr due 
to gravitational wave losses. Binaries with a ELM He-core WD are 
circled. Individual sources of interest are labeled: NLTT 11748 (the 
first eclipsing system ), CSS 41177 (the second eclipsing system; 
IParsons et all l20i"ll) , SDSS J0106-1000 and SPSS J0651+2844 
(two recently discov ered short period binaries; Kilic ct al. 120115 
IBrown et aT]|2011bl ). and KPD 1930+2752 (a binary containing 
a massive WD and a sdB star with a total mass e qual to the 
Chandrasekhar mass; IMaxte d. Marsh. lFNorthl 12000). The line 
at the Chandrasekhar mass M^h = 1.4 M© is where a merger 
cou ld produce a type la supernova in the tr aditional sc enario (but 
see van Kcrkwiik, Chang, & Justham 2010). Also sec Kilic ct al. 
(2012). 

stable ma ss transfer than that originally fou nd for cold 
He WDs (jMarsh. Nelemans. & SteeghsN2004l ). We per- 
form the full evolution calculations in §4, highlighting 
the new phases of H mass transfer, the special behavior 
near the period minimum, and the likelihood that the 
intruiging object HM Cancri is one of these systems. We 
close in §5 by discussing the implications for AM CVn 
evolution, and the remaining work needed to resolve the 
thermonuclear outcomes for the accreted matter. 

2. ELM STRUCTURE AND RESPONSE TO MASS LOSS 

The way in which the donor radius changes with mass 
loss differentiates the evolution of these ELM dono r bi- 
naries from the earlier work of iMarsh et al.l (]2004f ) for 
cold He WDs, where the WD always becomes larger 
as mass is removed. This is usually expressed as the 
logarithmic derivative of radius with respect to mass 
Cho = d\nRn c /d\nMn e which is « —1/3 for cold, de- 
generate matter. However, the thick outer layer of non- 
degenerate and stably-burning hydrogen in an ELM WD 
dramatically changes £hc- 

For this initial exploration , we u se the models 
of iSteinfadt, Bildsten. fc Arrasl (|2010al) that provide a 
range of hydrogen envelope masses M cnv for each to- 
tal mastQ, Mho- These represent a large range of ELM 



4 Some calculations such as Panei et all i 20071) and 
lAlthaus . Scrcnclli, & Bcnvc mxtol 1)2001 ) find that for masses 
greater than ^ 0.2 Mq, hydrogen will not burn stably but will 
instead undergo flashes. However, the precise boundary is not yet 
known and may be metallicity-dependent. For inst ance, the re- 
cently discovered binary SDSS J065133.33+284423.3 l|Brown et al.l 
1201 lbT l has an ELM WD with a high temperature and a radius 
indicative of burning even though it is somewhat more massive 
(0.25M©) than the expected burning limit, although for this 



a 0.04 



0.15/2.12 

0.15/2.51 

0.15/2.90 

0.15/3.30 

0.15/3.69 

0.15/4.09 

0.15/4.50 

0.15/4.91 

0.15/5.30 




of 1 Gyr, t h e exp ected envelope masses according to 
are 6 x 1CT 3 M (M He = 0.16 Mq) 
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Fig. 2. — Radii of M He = 0.15M© He WDs as a function of re- 
maining mass as the mass is stripped off. The models are labeled 
by different starting envelope masses, M e 

nv ) in units of 10 -3 Mq. 
We also plot the point where a mass equivalent to the envelope 
has been stripped (circles) and where the hydrogen mass fraction 
X = 0.01 (stars). Finally, the da shed line is Eggleton's zero- 
temperature mas s radius relation (jVcrbunt & Rappaporlj 119881 ; 
IMarsh ct al. 2004). Because of the diffusive tail, non-zero amounts 
of hydrogen extend well inward of the point where the stripped 
mass is equal to the envelope mass. 

WD models (0.125 M < A/ Hc < 0.20 M ) that ad- 
equately cover the possible envelope masses (M env = 
(1 — 5) x 10~ 3 M©) at the time of Roche lobe con- 
tact, or equivalently have ages of 1-13 Gyr; for an age 
of 1 Gyr, t h e exp < 
iPanei et alj ([20071 ) 

to 7 x 1O~ 4 M (M Hc = 0.25 M Q ), but the exact map- 
ping of M cnv to age is not known in detail. Instead, we 
verify the starting models by noting that they span the 
observed range in surface grav ity and effective tem pera- 
ture seen in ELM WDs (e.g., iVennes et al.l 1201 II ). with 
T cff = 8,000-20,000 K and log(g) =5-7. 

We need to find Cho as a function of the amount of 
mass that has been lost, AMn e , for each distinct initial 
model. These models have the temperature, pressure and 
abundance profile of a stably-burning pure H envelope 
in diffusive equilibrium with the He core. The result- 
ing transitional H/He layer has an abundance profile set 
by chemical equilibrium in t h e chan ging electric field, as 
described in ISteinfadt et al.l (|2010a|) . As we show later, 
the thermal time at the base of the H envelope is longer 
than the mass transfer timescale, so that the material re- 
sponds adiabatically to mass loss. To simulate mass loss, 
we simply go to the mass coordinate m = M to t ~ AMh c 
in the initial model and force that mass element to the 
surface by making its pressure artificially low (we went 
to P = 10 10 dyne cm" ; in contrast, the pressure at the 
H/He boundary is 10 17 dynecm~ 2 , while the pressure at 
the core is > 10 20 dynecm -2 ). 

As the underlying fluid elements now have lower pres- 
sures, we conserve their entropy by assuming that TP -7 
is constant, where the adiabatic index 7 was ~ 2/5 
for the region in question. The density is then com- 
puted from T and P f ollowing the equation-of-state (see 
ISteinfadt et al.ll2bl0al for details), and the composition 
of each mass element is assumed to be the same as in 

particular object the unknown contribution of tidal heating 
may contribute to its large size. We limit our exploration to 
M H e<0.20M©. 
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Fig. 3. — Response of the initially stably burning donor to mass 
stripping fjj e = dlnRn c /d\n Mjj c as a function of stripped mass. 
The initial models were for Mjjo = 0.15Mq He WDs with either 
a Menv = 0.003Mq or M cnv = O.OO4M H envelope, undergo- 
ing stable hydrogen burning. The dashed/dot-dashed lines show 
the results of this work using an adiabatic approximation to re- 
calculate the stellar structure, with the filled circles indicating the 
envelo pe mass. The soli d curves show the dynamical results of 
MESA {Paxton ct al. 2011) for a similar model (albeit with a lower 
envelope mass of M env = O.OO2IM0), using mass-stripping rates of 
10 _6 M Q yr _1 (blue), lO^Meyr -1 (green), and 10 -10 M Q yr -1 
(red) . Except for the transients at the beginning (which have to do 
with the outermost layers of the WD) all of the models generally 
agree for the first 10% of the WD. For the lowest accretion rate 
(not expected to be adiabatic) the model diverges after 0.015 M@. 
In all cases, fjje remains > 1 past twice the envelope mass. 



the original model. Therefore, we need only integrate 
the pressure and radius with respect to the mass coor- 
dinate. This integration is done by shooting from the 
core toward the surface and from the surface toward 
the core and meeting at a mass coordinate in the mid- 
dle (the X = Y point, where the hydrogen and helium 
mass fractions are equal). Selected results for 0.15 M Q 
ELM WDs are shown in Figure [21 where we plot the 
radius as a function of the remaining mass along with 
the zero-temperatu re (fully degenerate) model used by 
iMarsh et al.l (|2004p. Two things are obvious: (1) as the 
envelope is st ripped, the slope Cue > 1 so the donor 
shrinks (as in ID'Antona et al.l 120061 for stable burners) , 
and the slope is still greater than the —1/3 value of the 
zero-temperature model even once X < 0.01; and (2) 
the radii are always larger than the zero-temperature 
model, even when the envelope has been stripped, be- 
cause of the lower initial degeneracy. In our models the 
larger radii are consequences of both finite ages so that 
even just passi ve cooling would re sult in warmer/larger 
WDs (unl ike in IMarsh et al.ll2004D and stable H burning 
(unlike in Dcl ove et al.ll2007l) . The initial core tempera- 
tures at the onset of mass transfer range from 6.7 x 10 6 K 
for the 2.1 x 10~ 3 M Q envelope to 1.4 x 10 7 K for the 
5.3 x 10~ 3 M Q , with h igher temperatures resu lting in 
higher initial entropies (Dclove & Bildstcn 2003). 

In Figure [3] we compare our results for (ne with 
those computed using Module s for Experiment s in 
Stellar Astrophysics (MESA; iPaxton et al.ll201lD. For 
our starting points, the models of lSteinfadt et alj (|2010a) 
only covered limited values of Mh 6 and M cnv , as we used 
models that "evolved" from the starting models. Sim- 
ilarly, with MESA we had to use models that we could 



construct out of suitable starting conditions, subject to 
mass loss and diffusive equilibrium. Therefore the com- 
parison is not exact: the MESA model has a total hydro- 
gen mass of 2.1 x 10~ 3 M Q , radius of 0.048 Rq, and sur- 
face temperature of 7150 K, compared to (4.0 x 10 -3 M , 
0.062i? Q , 7160K) and (3.0xlO- 3 M Q , O.O48i? , 6500 K), 
for our models. We cannot simultaneously match the en- 
velope mass, radius, and surface temperature. Even so, 
the agreement for £hc is good, and the model profiles (the 
pressure and adiabatic index as functions of radius) agree 
well. The main apparent difference is the "hump" in £h c 
around a stripped mass of (2 — 3) x 10 -3 M Q visible in our 
models but not in the MESA models. This occurs slightly 
before the hydrogen-to-helium transition and near the 
onset of degeneracy in the original model, and given the 
different mass coordinates at which this happens in our 
models versus the MESA models, the discrepancy is not 
surprising. We see that overall, the shape of Ciie in this 
adiabatic limit is constant in terms of AMH e /M cnv . 

Whether or not adiabatic evolution is a valid assump- 
tion depends on the thermal timescale r t h at the base 
of the envelope, which is where the luminosity is fixed 
and most of the excess radius beyond the degenerate 
helium core is located. Mass-loss will be adiabatic if 
Menv/ 1 Ml < 7th- For the models of ISteinfadt et al.1 
(|2010al) . T t h is typically a few Myr, so for accretion that 
is faster than M cr i t = M onv /r t h ~ few x 10~ 9 M Q yr" 1 , 
our approximation is reasonable, and this regime encom- 
passes most of the evolution we discuss below. For lower 
mass-loss rates, we should consider thermal evolution 
during the mass stripping, although at these rates the 
amount of mass stripped is low enough that it does not 
significantly affect our results. 

3. MASS TRANSFER BASICS 

We start with a He WD donor with mass Mh c that 
is steadily losing mass onto an accretor (typically a CO 
WD) with mass M a . The total mass is M t = M a + 
Mne, while the binary orbital separation is a and the 



orbital angular momentum is J = M a Mn c {Ga/M t 



jl/2_ 

The mass transfer rate is fixed by the rate of orbital 
angular momentum loss, J via 



J M a , Mj 



He 



J M a M He 



a 
2a 



M t 
2M t 



(1) 



The accretor may have intermittent periods of unstable 
hydrogen burning (classical novae, or CNe) that eject 
material from the binary; therefore we cannot assume 
that mass transfer is conservative but must instead track 
the mass lost by the binary, M t . We write this as M t = 
/Mhc, so that / = 1 means the accretor, on average, 
keeps a constant mass (e.g., it ejects in each CN the 
amount of matter that has accreted), and f — means 
that the accretor keeps all the accreted mass. We do 
not recalculate the accretor 's radius as it gains mass, as 
this is a small effect for the massive accreting WD. If the 
CN were to excavate material from the accreting WD, 
then / > 1, although we do not explicitly consider this 
situation. We write the expression for J as 



J 

7 



Mi 



He- 



Mi 



He 



1 + (/-!) 



Mm 

M a 
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2M f 
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allowing a connection between the mass transfer rate and 
the loss of orbital angular momentum. 

We always assume that the donor's radius J?h b tracks 
that o f the Roche lobe, i?L, and use the simple IPaczynskil 
(fl967ft formula for R h , giving 



J 
J 



Ah 



He 



M He 



Chc 



+ (/"!) 



M Hc 



2M t 



He 



2 U j M a J 3M t 



(3) 



Here, £ rL = dlnTx/dln Miie is the derivative of the 
Roche lobe radius in units of a, tl = Hl/o. The 
Paczyhski approximation is co nsistent to with in about 
5% with the results using the lEggletorJ (|1983l ) formula 
for tl for the mass ratios considered here. 

At the onset of mass transfer, gravitational- wave losses 
set j to be (jLandau fe Lifshitz|[T975t) : 



Jgr 32 G 3 M t M Hc M a 
J 5 c 5 a 4 



(4) 



This remains true for wide systems, as the material 
transfers through an accretion disk. However, if there 
is not enough room for a disk (based on Eqn. 6 of 
lNelemans~e t al. 20( 53) then material will impact dire ctly 
onto the accretor (|WebbinkllT98l iMarsh et al.ll2004t ). In 
this case, the angular momentum of the accreted ma- 
terial is lost to the accretor. This angular momentum 
loss subtracts an additional \Avi(l + M^ c /M a ) from the 
quantity in the brackets in Eqn. ([3]), where ry t is the 
equivalent radius of the material orbiting the accretor 
(in units of a) as given bv lVerbunt fc Rappaportl (|1988l 
correcting for their inverted definition of the mass ra- 
tio). In reality, all of the J may not be lost but some 
may be transferred back to the system via t idal coupling 
(as in IMarsh et alll^OOl iFuTler fc Lai|[20lJ) . but for now 
we explore the limit where this does not occur since the 
tidal timescales are largely unconstrained for the accret- 
ing WD. 

3.1. Burning of the Accreted Material 

The importance of knowing the value of / is high- 
lighted by Equation ([2]). When / = (i.e., the trans- 
ferred material stays on the accretor), the mass ratio 
Mn c /M a decreases even more quickly, reducing dynami- 
cal stability (§ 13. 2|) and increasing |Mhb|- On the other 
hand, when material is lost from the system (i.e., / = 1), 
then Mn e /M a changes more slowly with a corresponding 
reduction in |-Mne|- It is the ability for the accretor to 
burn the accreting material at the supplied rate that de- 
termines /. Stable burning (i.e., burning the material 
at the accreted rate in a thermally stable manner) im- 
plies / = 0, whereas unstable burning results in mass 
losing events (e.g., classical novae) that drive / — > 1. 
Since no detailed calculations are available that cover ac- 
creted material with our particular abundance mixture 
(almost pure H going to nearly pure He); we base our cur- 
rent c a lculations on IPapaloizou. Pringle. fc MacDonaldl 
(fl98l^ : lNomoto et all (120071) :IShen fc Bildstenl (l200lT fo"r 
H burning and llben fc Tutukovl (|1989ft "for" pure He burn- 
ing. 
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Fig. 4. — Evolution of compact WD binaries involving an ELM 
He WD and a 0.8 Mq CO WD. The masses of the He WDs vary 
from 0.125 Mq to 0.2 Mq, and the envelope masses (in units of 
10 — 3 Mq) are also listed. We plot the expected mass-transfer rate 
M from Roche-lobe overflow as a function of orbital period after 
the systems are brought into contact by GW emission. The sys- 
tems begin on the lower tracks, move toward shorter orbital periods 
and increasing M, reach a period minimum, and then leave along 
the upper tracks. The numbers show the ages of the systems (after 
contact) in Myr. The thick regions are where the accretion stream 
hits the CO WD directly, instead of forming an accretion disk. The 
tracks are additionally labeled with symbols that show the hydro- 
gen fraction X of the accreted material, changing from 99% down 
to 0.1%. The black das hed line is Eqn. El The r egions where stable 
H UNomoto e t al. 2007) or He l llben & Tutu kov 1989) burning are 
possible are shown by the hatched regions, although we use the 
green horizontal line for the lower H-burning limit to account for 
reduced metallicity. Finally, the red dashed line shows where ac- 
cr etion is unstabl e to dwarf novae (dNe) for pure hydrogen, based 
on lShafterl l[T99l) . 

The transferred material at the onset of mass transfer 
is certainly pure hydrogen. If that material does not gain 
any carbon or oxygen contamination from the accreting 
WD, then it can on ly burn via the pp-cycl e. Though 
briefly considered by IPapaloizou et ail (fl982) . no recent 
calculations have been done to explore the thermal sta - 
bility of such a pure H layer. iShen fc Bildstenl ([20071 ) 
performed stability calculations for very low metallic- 
ity, going down to 10~ 4 Solar metallicity. They found 
that the burning material was more stable, reducing the 
minimum accretion rate where stable burning can oc- 
cur by over a fac t or of three from the limits derived 
by iNomoto et all (|2007l Eq. 5) for Solar metallicity 
(shown by the the red banded region in Fig. 2]). We 

therefore assume stable hydrogen burning when M > 
10" 7 (M a /M Q - 0.5357) Moyr" 1 (the green horizontal 
line in Fig. 2]), and assume that all the accreted mass re- 
mains on the accretor (/ = 0). For rates lower than this, 
we assume that the H burning is unstable, leading to 
Roche lobe overflowing nuclear flashes and / = 1. How- 
ever, very little mass is transferred at these low rates, 
leading to a small impact on the outcome. 

The accretion rates needed for stable pure He burn- 
ing are much higher than those for hydrogen. Just 
as in H burning, the stability arises due to the in- 
creasing importance of radiation pressure and shell 
thickening, both occurring as t he burning luminosity 
approaches the Eddington limit (|Tutukov fc Yungelsonl 
ll996tlShen fc Bildstenll2007f) . Most thoroughly discussed 
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in llben fe Tutukovl (|1989| ). this rate depends on the WD 
core mass, and, for the evolution of these binaries, al- 
lows for stable He burning on ly near the period minimum 
(jTutukov fe Yungelsonl[l996| ). This results in a narrow 
range of rates (typically only a factor of ~ 2, see the 
hatched region in Fig.[4j for which a WD can stably burn 
helium. This demands two comparisons for our calcula- 
tions. First, a tracking of when the accretion rate is in 
the stable regime and for that duration, a confirmation 
that enough mass is transferred to build the steady state 
burning model. Second, when M reaches the upper limit 
of the stable burning regime (meaning that the radius of 
the accretor exceeds its Roche radius), we must account 
for the excess transferred matter that cannot be stably 
burned by the accretor. For t his calculation, we s imply 
assume (as is commonly done: lNomoto et al.l 120071 ) that 
that matter is lost by the binary. Hence, even in the 
stable regime we will sometimes set / to a value that is 
different than zero0 We do this so that |MHe| does not 
exceed the limit for stable He burning (essentially the 
Eddington limit). When these conditions occur in our 
evolutionary scenarios, its almost always the case that 
X R3 0. Eventually, as the orbit widens, |Mh | decreases 
below the lower limit for stable He burning, and we set 
/ = 1 again. 

3.2. Enhanced Stability of Mass Transfer from an ELM 

Mass transfer is dynamically stable when the term in 
the brackets in Eqn. © is positive for / = 0. If it is 
negative, the stellar radius will be larger than the Roche 
radius after mass transfer, leading to a dynamical insta- 
bility and merger. The form of this condition depends 
on whether the angular momentum of the material which 
leaves the donor can get back into the orbital angular mo- 
mentum (as assumed in disk accretion) or becomes "lost" 
in the accretor during direct impact. Accretion via a disk 
is stable when 

Mac 5 Chc ,r\ 

MT < 6 + — ' (5) 

For a cold WD, i?Hc ^ A/ Hc ' > leading to ^He = —1/3. 
This would imply stability if Mn c /M a < 2/3, the tra- 
ditional criterion for stability. The additional angular 
momentum loss allowed by direct-impact changes the cri- 
terion to: 



Afac 5 
M a 6 



CHc 

2 




(6) 



This substantially reduces t he allowed mass ra tios for ini- 
tially stable mass transfer ([Marsh et al. 2004) especially 
for cold WDs, or those with a finite entropy ([Delove et alJ 
2007) but no burning. 

However, as we show in Figure [3l He WDs with thick, 
massive hydrogen layers h ave Ch s ^> 1 for str ipped 
masses AMho/Mhc < 10%. D'A ntona et al.1 ((2006T) were 
the first to explore the possibly different evolutions this 

5 In contrast to the majority of our numerical integrations where 
we use an adaptive step size 4 th /5 th -order Runge-Kutta integrator, 
updating the value of / is done so as to keep the rate of accreted 
material for the last iteration in the stable regime — no higher 
order derivative is used — which leads to the apparent numerical 
noise in Figure [4] This does not affect the outcome. 



allows for. Just from Eqn. ([3]), Mn e < M a mass transfer 
will be dynamically stable as long as there is an accre- 
tion disk. Even without a disk (using Eqn. [6]), the ad- 
ditional term + Mft c /M a )rh is roughly ps 0.5, so as 
long as Chc is positive we still have dynamical stability. 
Hence, for almost all the mass ratios considered here, dy- 
namical instability is not possible for disk-fed accretion 
and could only happen for direct-impact accretion when 
Mn e /M a > 1/4 (assuming a typical £hc = —1/3 once 
the hydrogen envelope has been fully stripped). Tidal 
interactions could work to stabilize these f ew unstable 
mass transfer cases (as in lMarsh et al.ll2004D . but that is 
beyond the scope of our work. 

4. EVOLUTION INCLUDING MASS TRANSFER FROM AN 
ELM WD 

Just as in the conventional double WD binaries, the 
lower mass ELM WD will be driven into contact by 
the loss of angular momentum from gravitational waves. 
For the more massive companion, we presume cold CO 
WDs with masses 0.7-1.0 M Q (we do not discuss masses 
< 0.7 M^, as those ar e not covered by the calculations 
of lNomoto et al.l 120071 ). We neglect any tidal heating of 
the ELM WD prior to contact, and simply assume, as in 
the previous section, that it has a thick H envelope that 
is stably burning. We start by discussing the long pe- 
riod of stable mass transfer of the overlying H shell, then 
discuss the behavior near the period minimum and close 
with a brief discussion of the outgoing He mass transfer 
phase; which is similar to that previously discussed in 
the literature. 

4.1. Mass Transfer Phase up to the Period Minimum 

The larger stellar radius leads to the onset of mass 
transfer at orbital periods of ~ 10 minutes, much longer 
than expected for a thin-H-shell cooling WD. Since 
Chc ^> 1, the donor shrinks under mass loss, so the mass 
transfer is stable and secularly driven at the rate implied 
by equation ([3]) when angular momentum loss is fixed 
by gravitational radiation (see Eqn. |4j. We explore a 
range of ELM WDs with masses M Rc /M Q =0.125-0.20 
and envelopes of M cnv = 1.5 x 10~ 3 to 5.1 x 10~ 3 M Q , 
although a single donor mass only has part of the range. 
Not all of these are physically realistic: some higher-mass 
WDs may not have stable burning, or some combina- 
tions of envelope and total mass m ay not be reachab le in 
a Hubble time (see di scussions in lPanei et alJ [20071 and 
ISteinfadt etllll2?il0ah . 

We st art our calculation wh en mass transfer begins 
(also see iD'Antona et ah! 120061 ). typically at orbital pe- 
riods of 6-10 minutes. Figure [5] shows the resulting evo- 
lution up to the period minimum for two cases with 
Mne = 0.15 M Q and 0.2 M®. For these cases, the ac- 
cretion rate, M , is in the range where our adiabatic ap- 
proximation is valid0 The increase in \M\ as the or- 
bital period shrinks is evident. At the earliest part of 
the onset of accretion, the orbital P (second panel) is 
nearly identical to that calculated (dot-dashed lines in 
the second panel) from the loss of angular momentum 

6 Other cases we show in a Figure [4] that come into contact at 
much wider orbital periods violate this inequality, but in those 
cases the amount of mass transferred at this low rate is small 
enough as to not qualitatively change the results. 
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from two orbiting point masses. Namely, if measured 
at this stage, the system would appear to have an or- 
bital period change co nsistent with gravitational inspiral 
(jD'Antona et al.ll2006l also see Fig. [6]) . During this initial 
evolution we have £ne 3> Cr L > an d we find from Eqns. (|3j) 
and (jU that 



|M He 



64Af£ c M a /256tt 8 G 5 



5c 5 Chc V MtP 8 



1/3 



(7) 



highlighting the need to calculate the considerable evo- 
lution in Chc during this period (Fig. [5j to find Mh c - 

This evolution is a new mode of H mass transfer onto 
a WD that has not been previously explored over a wide 
range of initial masses with semi-analytic models. The 
long (~ Gyr) life of the ELM WD prior to contact has 
allowed for complete diffusive equilibrium so that the 
transferred material is pure hydrogen at the start, with 
an increasing amount of helium at later times as the 
stripping reveals the underlying He WD core (5th panel 
in Figure [5}. The amount of mass that has been stripped 
(third panel down) is simply what's needed to keep the 
stellar radius equal to the Roche radius as the donor 
moves inward. These are often but not always disk accre- 
tors but compact enough that the thermal disk instabil- 
ity that gives rise to dwarf novae is suppressed (systems 
below the dashed line on Figures 0] and [7] would undergo 
dwarf novae instabilities). 

The thermonuclear outcomes during this phase remain 
uncertain until a time dependent accretion calculation 
has been performed with the changing M and X of Fig- 
ure [5] However, as we discussed in § 13. H if H mass trans- 
fer leads to stable burning, the system's luminosity would 
increase to L = Q nuc M. 

4.2. Near the Period Minimum and The Ultra-Compact 
Binary HM Cnc 

After the initial evolution the degenerate portions of 
the donor come to the surface and a minimum period is 
reached. This minimum is near where \M\ is maximum, 

where direct impact often starts (forcing us to adjust J 
accordingly), and where the mass transfer can in princi- 
ple become dynamically unstable. For less massive ac- 
cretors direct-impact accretion will always happen. The 
mass-transfer rate is shown in Figures 0] and [7] for a sin- 
gle accretor mass and a range of donor masses, while 
Figure [6] shows the orbital period derivative. The mass- 
transfer factor / is calculated according to § 13.11 

The interacting binary star HM Cancri (HM Cnc 
or RXJ0806 . 3+152 7 ) has an orbital period of 5 . 4 min 
(llsrael et all 119991 : iRamsav. Hakala. fc Cropper] 120021: 
llsrael et all I2002D . Spectroscopy (|Roelofs et all 120101) 
supports a model with two interacting white dwarfs con- 
sistent with masses of 0.27 M Q and 0.55 Mq but with the 
actual values unconstraine d. X-ray observati ons show a 
decreasing orbital period (jStrohmaverl [2005) at a rate 
consistent with gravitational-wave emission, as shown in 
Figure [6] (where we highlight it's value of P and P). 
While the X-ray luminosity was puzzli ngly low for ac- 
cretio n from a degenerate companion, iD'Antona et al.l 
(2006) showed that the implied transfer rate was much 
more consistent with an ELM WD donor. We find the 
same conclusion, which is that HM Cnc is consistent 




0.1 0.01 0.001 

Time Before Period Minimum (Myr) 

Fig. 5. — Evolution of compact WD binaries involving an ELM 
He WD and a M a = O.8M CO WD. The masses of the He 
WDs arc A/ Hc = 0.15M© (blue solid line) and M He = 0.2M Q 
(red dashed line), with envelope masses Of Menv = 2.1 X 1O" 3 M 
and A/ onv = 1.8 X 10~ 3 Af , respectively. We plot (top to bottom) 
the orbital period, the orbital period derivative, the transferred 
mass, the accretion rate, the hydrogen mass fraction, and (jie- 
The period of direct-impact accretion for the Ainc = O.2A/ WD 
is indicated by the thick dashed line. In the P panel, the red and 
blue dot-dashed lines are where the loss of angular momentum is 
entirely given by Eqn. J4f. 




4 5 6 7 8 9 10 

Orbital Period (min) 

Fig. 6. — Period derivative vs. orbital period for a compact WD 
binaries involving an ELM He WD and a 0.8 M CO WD. The 
different donor models are as in Figure [4] We also the expected 
P oc — p -5 / 3 behavior along the incoming (lower) track (Eqn. [Tol l 
as expected from gravitational inspiral. We also show the observed 
properties of HM Cnc. 



with accretion from an ELM WD that is still inspiral- 
ing and transferring H-rich material. Such a system, 
along with oth er short-period binar ies like V407 Vul (RX 
J1914.4+2456; lCropper et alj|1998ft . are extremely valu- 
able as probes of the eventual fates of systems like those 
in Figure [T] 

4.3. Impact on AM CVn formation 
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Differentiating between the different possible outcomes 
based on the progenitors characteristics is a major goal 
of work like that described here. In particular, we seek 
to separate the systems that will remain separate from 
those that will merge. Systems that remain stable will 
move to longer orbital periods with declining accretion 
rates and become AM CVn binaries: systems that have 
a prolonged period of He mass transfer driven by grav- 
itational wave emission. In contrast, those that become 
unstable will shred the He WD into a rapi dly rotating hot 
disk t hat likely becomes a n R CrB star (|Saio fc Jefferyl 
12002ft . IMarsh et all (|2004ft did this in detail by consider- 
ing the stability of mass transfer, including the unknown 
effects of tidal coupling. 

Based on the qualitative effects described in § 13.21 we 
can examine here whether more systems will remain sta- 
ble. This will directly in fluence population calculations 
(e.g.. lBrown et aDl2011af ). which find only a small frac- 
tion of the total AM CVn population could come from 
ELM WD binaries, although based on the limited param- 
eter space that we explore we do not make population 
predictions in this work. 

In Figure|8]we differentiate between those systems that 
has disk-fed accretion and those that undergo direct im- 
pact. For all of the systems that we study dynamical 
stability (Eqns. [5] or is maintained. We also see that 
the increased size and decreased degeneracy of the donor 
which begins mass transfer at larger separations means 
that more space is available for accretion disks. The ex- 
act threshold depends on both the donor's core and en- 
velope masses, although in general it increases by about 
0. 03 Mp, or so compar ed to the threshold for a cold donor 
in IMarsh etal] (|2004D . 

While none of the systems we study is dynamically 
unstable, this is largely because we do not include the 
less massive accretors (0.4-0.6 M©) which appear com- 
mon (jKilic et al.ll2012f ) and which would be closer to the 
stability threshold. Extending our calculation to lower 
accretor masses with the addition of improved stable- 
burning thresholds will help investigate dynamical sta- 
bility more closely, but the qualitative behavior of an in- 
creased range of systems maintaining accretion disks and 
remaining stable will hold. Current estimates, based on 
rather small samples of ELM WDs and AM CVns, sug- 
gest that o nly a small fraction of AM CVns come from 
ELM WDs (|Brown et al.ll2011af ). but signif icant imp rove- 
ment s in the AM CVn population (e.g., iLevitan et al.l 
120111 ) and evolutionary calculations may change this con- 
clusion. It is clear that this increased phase space for 
stable accretion outcomes demonstrated here will lead to 
a higher yield of AM CVn binaries from double WDs. 



4.4. Stable Outgoing Helium Transfer 

The final behavior is familiar and the change in J is 
determined not by the change in a but by the mass trans- 
ferred. The mass-radius relation of the donor determines 
the evolution. In this range, £jje is roughly constant (the 
classical value is —1/3, and in our model is it close to 
-1/4). With that, we find 



Po 



Mi 



He 



M 



Co 



M, 



-1/3 
He 



3/2 



(8) 



where i?Hc = Ro{Ma e /Mop a , and Co is the final constant 
value of Cho- With this, we ignore the terms in the brack- 
ets in Eqn. ([2} that contain Mue since they are less than 
the other terms (5/6 + (jje/2), and find: 



|M Hc 



7776 M M a B^ ( '4tt 14 G 



i4r<2\ l / 3 



V AftP" 



(9) 



in the limit that Co ~ -1/3. With Co = -1/4, the ex- 
ponent on P changes from —14/3 = 4.67 to —104/21 = 
4.95. 

Based on these, we can determine P along the incoming 
and outgoing branches: 



Pn = 



-384Af Hc /47r 8 G 5 M a 2x 1/3 



Pn 



5c 5 V p 
1728 {2tt 22 G 7 R 9 M 3 M*\ 1/6 



5c 5 



P 



10 



(10) 



This then allows us to determine the relative populations 
of objects along the incoming and outgoing branches. 
This is just the ratio of the P's, since the time spent 
(oc 1/P) is related to the population: 



out fx 



p, 



Pn, 



^8G PM He 



3/2 ' 



M o 9nR 3 



(11) 



Based on our model, R = 0.034P© at M = 0.05M©, so 
^out/in ~ 4(P/10min)(M He /0.15M©) where the mass 
of the He WD is the original mass along the incoming 
track. This shows that if every AM CVn had an ELM 
WD as the originating double WD binary, then there is 
one inspiraling H mass transferring system for every 4 
AM CVns at an orbital period of 10 minutes. Note that 
this is not true in the period/mass range where direct- 
impact accretion occurs (as presumably for HM Cnc), 

as that dramatically increases |P ou t| and so reduces the 
outgoing (AM CVn) population. In any case, finding 
t hese predicted sy s tems remains a challenge. 

ID'Antona et~aT1 (2006) found it unlikely that V407 Vul 
had the same origin as HM Cnc, as the lifetimes on the 
incoming versus outgoing branches would suggest that 
its orbital period should be increasing rather than de- 
creasing (|Strohmaverl 120021 find P = -2.6 x lO^ss -1 
at an orbital period of 9.5 min). We do not find quite 
such a strong preference, with A/" ou t/in ~ 2-4 depending 
on the accretor mass (c ompared to a ratio of about 5 
in ID'Antona et aT] 120061 ). Still, the small orbital period 
derivative of V407 Vul is somewhat difficult to explain, 
being a factor of 2 smaller than the lowest value we find 
in Figure [6] Considering lower accretor masses may help 
resolve this. 

5. CONCLUSIONS 

We have shown that the unique properties of ELM 
WDs — large, nondegenerate H-rich shells supported by 
stable H burning — lead to some new phenomena when 
mass transfer initiates in double- WD binaries. There 
is a prolonged period of H-rich mass transfer at a low 
rate during inspiral, with HM Cnc potentially being the 
prototype. The change in the mass-radius relation for 
the donor creates an intrinsically more stable binary 
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birthrate, potentially alleviating the apparent paucity of 
progenitor systems. Those AM CVns which emerged 
from this progenitor scenario will also have a larger He 
core radius than expected from an initially cold WD, 
thereby exhibiting a higher accretion rate at a fixed or- 
bital period than from a cold WD (jDelove fc Bildstenl 
l200l . 

Prior to the onset of mass transfer, the ELM WDs had 
Gyrs to undergo diffusive settling and substantial burn- 
ing of hydrogen. That clearly allows for the complete 
sedimentation of the heavier elements from the outer- 
most layers of the WD. Hence, the mass transferred will 
vary from nearly pure H, to nearly pure He. As we dis- 
cussed, much work remains to more carefully calculate 



#D ^rjmp^uglgr autaojngft. fgfifi3eihjg: tp%3^ei e f^f 
more thermally stable, then these systems may become 
more observationally detectable due to the higher lumi- 
nosities. It is also interesting to note the pronounced ab- 
sence of heavy elem ents in the x-ray spectra of HM Cnc 
(|Strohmav cr 2008), also pointing to an ELM origin that 
lived a long time prior to mass transfer initiation. If ther- 
mally unsta ble, then the accumu lated mass could ignite 
explosively (jBildsten et al.l 12007). potentially contribut- 
ing to the increasing number of low-luminosity " super- 
novae" observed locally (e.g.. iKasliwal et aT1l2010l ). 

We thank Bill Paxton for advice on running MESA, 
and thank Bill Paxton and Ken Shen for helpful dis- 
cussions. This research has been supported by the Na- 
tional Science Foundation under grants PHY 11-25915 
and AST 11-09174. 
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